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Abstract — We develop a two-dimensional (2D) Superconduct- 
ing Quantum Interference Filter (SQIF) array based on recently 
introduced high-linearity tri-junction bi-SQUIDs. Our bi-SQUID 
SQIF array design is based on a tight integration of individual bi- 
SQUID cells sharing inductances with adjacent cells. We provide 
extensive computer simulations, analysis and experimental mea- 
surements, in which we explore the phase dynamics and linearity 
of the array voltage response. The non-uniformity in inductances 
of the bi-SQUIDs produces a pronounced zero-field single anti- 
peak in the voltage response. The anti-peak linearity and size can 
be optimized by varying the critical current of the additional 
junction of each bi-SQUID. The layout implementation of the 
tight 2D array integration leads to a distinct geometrical diamond 
shape formed by the merged dual bi-SQUID cells. Different size 
2D arrays are fabricated using standard HYPRES niobium 4.5 
kA/cm 2 fabrication process. The measured linearity, power gain, 
and noise properties will be analyzed for different array sizes 
and the results will be compared with circuit simulations. We 
will discuss a design approach for the electrically small magnetic 
field antenna and low-noise amplifiers with high bandwidth based 
on these 2D bi-SQUID SQIF arrays. The results from this work 
will be used to design chips densely and completely covered in 
bi-SQUIDs that has optimized parameters such as linearity and 
power gain. 

Index Terms — SQUID, numerical simulations, small antenna, 
low noise amplifier, flux noise, high sensitivity 



I. Introduction 

THE quest to increase the linearity of SQUID and SQIF 
arrays was boosted by the introduction of the bi-SQUID 
- a tri-junction dc SQUID with a linear voltage response (TJ- 
(3). A non-linear inductance of the additional third junction 
provided the desired linearizing effect. These bi-SQUIDs are 
now being used in uniform and non-uniform (SQIF) arrays 
in place of conventional dc SQUIDs with a goal of achieving 

Manuscript received October 9, 2012. This work is supported in part 
by the Tactical SIGINT Technology Program N66001-08-D-0154, SPAWAR 
SBIR contracts N00039-08-C-0024, N66001-09-R-0073, ONR (Code 30), 
the SPAWAR internal research funding (S&T) program and HYPRES IR&D 
Program. 

S. Berggren and A. Palacios are with San Diego State University, San 
Diego, CA 92182 USA (e-mail: susan_berggren@yahoo.com). 

G. Prokopenko and O. A. Mukhanov are with HYPRES Inc., Elms- 
ford, NY 10523 USA (phone: 914-592-1190; fax: 914-347-2239; e-mail: 
georgy @ hypres .com, mukhanov @ hypres .com) . 

P. Longhini, A. Leese de Escobar, B. J. Taylor, M. C. de Andrade, E. Wong 
and V. In are with SPAWAR SSC Pacific, San Diego, CA 92152 USA. 

M. Nisenoff is with M. Nisenoff Associates, Minneapolis, MN 55403 USA. 

R. L. Fagaly is with Quasar Federal Systems, San Diego, CA 92121 USA. 

E. Cho and T. Wong are with the University of California, San Diego, La 
Jolla, CA 92093 



higher linearity. To date, the most of the design efforts using 
both SQUIDs and bi-SQUIDs went into the optimization of 
one-dimensional (ID) serial or parallel arrays and their com- 
binations. This is due to the higher complexity analysis and 
modeling required for 2D arrays to account mutual coupling of 
neighboring cells and complex current distribution in arrays. 

Arrays of DC bi-SQUIDs with size varying from loop to 
loop in an unconventional geometric structure are known to 
exhibit a magnetic flux dependent average voltage response 
< V(<& e /<&o) >, where <£ e denotes the external magnetic flux 
and $o is the magnetic flux quantum, that has a pronounced 
single peak with a large voltage swing at zero magnetic 
flux. The potential high dynamic range and linearity of the 
"anti-peak" voltage response render the array a promising 
candidate for multiple applications including a detector of 
absolute strength of external magnetic fields, wide-band low 
noise amplifiers and magnetic antennas. These arrays are also 
commonly known as Superconducting Quantum Interference 
Filters (SQIFs). Since it was theoretically proposed El, (5j 
and experimentally demonstrated for the first time [6]-| 8 ] the 
SQIF concept is investigated and exploited by a continuously 
growing number of groups with respect to its basic properties 
|9)-|[T4| and its suitability in various fields of application 
like magnetometry p5|-pQ|, low noise amplifiers of radio 
frequency (RF) and digital signals (2TJ-||29) and electrically 
small antennas (30|, (3T| . The unique noise properties, wide 
bandwidth and linearity of certain configurations of the SQIF 
array makes it especially attractive for an integration with su- 
perconducting wideband Digital-RF receivers being introduced 
for a variety of military and civilian applications |32|-p7). 
These receivers are based on highly linear superconducting 
analog-to-digital converters [38]-[40] which require equally 
high linearity of analog RF electronics. 

In this paper, we present a new design for the 2D array 
- a tightly coupled 2D network of bi-SQUID cells, in which 
junctions and inductances are shared between adjacent cells. 
In order to achieve this, we alter the shape (or inductance 
selection) of the bi-SQUID cells and merge two cells together 
to form a dual bi-SQUID - a diamond- shaped bi-SQUID. 
A diamond shape layout of these cells makes 2D array 
integration easier and allows for the bi-SQUIDs to be more 
densely packed. We demonstrate numerically and experimen- 
tally that these 2D bi-SQUID SQIF arrays can be optimized 
to produce a linear anti-peak at the zero magnetic flux. We 
examine in great detail the voltage response as a function of 
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Fig. 1. Schematic diagram of a single DC SQUID (left) and a single DC pig. 2. Average voltage response of a single DC SQUID compared with 

bi-SQUID (right) where (i b , h, i 2 , 13, U, 15) are the normalized currents, that of a sing i e DC bi-SQUID. The computer simulations of the system of 

Opi, ^2, ^3) represent the phases across the Josephson junctions, L/2 is the equations in Eq. {TJ and Eq. {2} were performed with i b = 2.0, L = 1.0 and 

parameter related to the inductance values and x e is the point in the equations j c3 = 1.0. 
where the external fields are included. 



controlled parameters, including: inductive coupling between 
loops, number of loops, bias current, and distribution of loop 
areas. The results from the numerical analysis are validated 
against experimental results from fabricated designs. 

II. Background: The single SQUID and bi-SQUID 

A single DC SQUID has two Josephson junctions arranged 
in parallel, connected with superconducting material, see the 
schematic diagram on the left in Fig. [T] The equations are 
derived using Kirchhoff's current law [41], which is the 
principle of conservation of electric charge and implies that at 
any junction in an electrical circuit the sum of currents flowing 
into that node is equal to the sum of currents flowing out of that 
node, along with a resistively shunted junction (RSJ) model of 
the over-damped Josephson junction. The Josephson junctions 
are assumed to be symmetric, in particular having identical 
critical currents I c \ and i C 2 that we use to normalize all the 
other currents in our consideration. The system of equations 
that models a single SQUID dynamics (42) is 



(fl = — - —{<pi ~(f2- <Pe) ~ Sm^i 

^2 = ^ + -^Oi - ^2 - (fe) ~ sin (^ 2 , (1) 

where (fi and Lp 2 are the phases across each of the Josephson 

junctions and the dots denote the time differentiation with 

2eI c R N 



1 



normalized time r = oo r t, where t is time and uo r 



h 



The parameter Rn in u c is the normal state resistance of 
the Josephson junctions, I c is the critical current of the 
Josephson junctions, e is the charge of an electron, and h 



is the reduced Planck constant. % 



h 

Ir 



is the normalized 

He 



bias current, where R is the bias current. L = 2tt—- is the 

$0 

B e 

normalized inductance and <z? e = 2nax e , where x e = — - is 

^ $0 
the normalized external magnetic flux per unit area and a is 

the bi-SQUID area. We use the approximate assumption that 
a = L. $0 = — « 2.07 x 10 -15 tesla meter squared is the 
flux quantum, wfiere h is Plank's constant and 2e is the charge 
on the Cooper pair. 

A DC bi-SQUID, which is a SQUID with an additional 
Josephson junction bisecting the superconducting loop, was 



introduced as an alternative to traditional SQUIDs and have 
shown superior linearity in the average voltage response anti- 
peak feature. The phase equations for the single DC bi-SQUID 
schematic on the right in Fig. [T] are derived in a similar way 
to those of the single DC SQUID. In this case there is a third 
junction J3 that is related to the first and second junctions 
through the phases: ipi + (fs = (f2- Using this relationship 
all the terms that include (fs can be replaced with ^2 — ^1, 
thereby reducing the number of phase equations needed to 
model the system from three to two as shown in (T). The 
governing equations for a single bi-SQUID are 
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sin (f2, 
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and I c s is the is the critical current of the 

1 c 

third junction. All other parameters are defined as the single 
DC SQUID. 

The single DC SQUID and DC bi-SQUID are simulated 
by integrating the systems of equations that model the system 
dynamics, Eq. ([T]) and Eq. ^ respectively. After the phases tpi 
and (f2 have been determined the derivatives (pi and (^2 are 

evaluated and the time-dependant voltage V(t) = 

is calculated. The average voltage, < V >, of a SQUID (or 
bi-SQUID) at a point in x e is the mean value of the voltage 
over time, presented in Fig. [2] The average voltage response of 
the bi-SQUID, with the proper selection of parameters, has a 
more linear average voltage response than the conventional DC 
SQUIDs. The higher linearity of the average voltage response 
increases the utility of the device as a linear amplifier. 

III. Diamond Shaped bi-SQUID 

The repeating pattern in the 2D array is a diamond-shape 
created by two bi-SQUIDs. Fig. [3] depicts a circuit of a single 
diamond structure of bi-SQUIDs where (ib, h,—Mi) represent 
the normalized currents, ((fi...(pe) are the phases across the 
Josephson junctions and (Li, L2 a ,—,^6a> ^2bv-^6fe) represent 
the normalized inductances. To assist with the design of the 
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Fig. 3. Single diamond structure of bi-SQUIDs circuit representation 
where ii,...,in) are the normalized currents, (<£i...<£6) represent the 
phases across the Josephson junctions, (L±, L2 a ,--,LQ a , L2b,---,LQb) are the 
parameters related to the inductance values and x e is the point in the equations 
where the external fields are included. 



array, equations that model the system are derived to be used 
in computer simulations test the effects of different parameter 
values. Six governing equations are needed to simulate the 
average voltage response for a single diamond structure of 
bi-SQUIDs, which contains six Josephson junctions and nine 
inductors. 

Kirchhoff 's current law results in the following relations for 
the currents and phases in the diamond- shaped dual bi-SQUID 



ib = n H- 
h + h = 
lg + h = 
hi i6 



^2 
^9 

hi 
h h 
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where are the phases across of the Josephson junctions, 
i = 1, 6. I c i are the critical currents of each of the junctions, 
z 1, ...,6, that satisfy I cl I c2 I c a h- hz 

Ica/Ic is me normalized critical current of the third junction 
and i C Q = I c q/I c is the normalized critical current of the sixth 
junction. The dot denotes time differentiation with normalized 



time r = u c t. Applying the mesh rule for the phase of the 
macroscopic wave function of the Cooper pairs to the top half 
of the diamond gives 

ipi + Lii 7 + L 2a i 8 = 27rx e ai + (p 2 + £26^9, 

where (p e — 2iraix e , where x e is the normalized external 
magnetic flux and we use the approximate assumption that 

a\ = Li+L 2a +L 2b . Combining ij = ig—ho and i 7 = hi—ig 

such that h = — — — + — — , and then substituting the 

current relations for i 7 , ho, i>9, and ig yields 



Li T V Li . Li . r . 

— + ^2a In ^"U + ^-^5 + ^1^6 =<f2 — <Pl 



L 2b i 2 + 27rx e ai + (Li + L 2a + £25) ^3- (4) 



In order to get the first of the six equations that is needed to 
describe the dynamics of the diamond shape, the bias current 
relation i 2 = i b — h and Josephson junction relations from 
Eq. ([3]) are substituted into Eq. ^ to become 



L 2b h 



Li 2 ((fi - 3 ) + y (20 6 + 05 - 04) = ^ y 

+ <p 2 - (fx + 27rx e ai + L12 (i C 3 sin <p 3 - sin 

+ y ( sin ^4 - sin (p 5 - 2i c6 sin cp 6 ) , (5) 

where L12 = Li-\-L 2a -\-L 2b . To solve for the second of the six 
equations that governs the dynamics of the single diamond the 
current relations for the bias current h = h~h and Josephson 
junctions from Eq. Q are substituted into Eq. ^ to give 



Li 



(20 6 + 05 - 04) - £l2(02 + 03) 



Li 



L 2a h 



+ (f 2 - (fi + 2irx e ai + L12 (z c3 sin <p 3 + sin <p 2 ) 
+ (sin y? 4 - sin <p 5 - 2i c6 sin <p 6 ) . (6) 

The current-phase relation around the upper loop in the single 
diamond is 

<Pl + (f3 + ^3^3 = <£>2, 

where L3 = + L35. Substituting the Josephson junction 
relations from Eq. and reorganizing yields the third of the 
six equations for the dynamics of the single diamond. 



^303 = ^2 - ipi - ^3 - L 3 i c3 sin(^ 3 . 



(7) 



A similar procedure is used to determine the equations for 
the three junctions located in the bottom half of the diamond 
structure and can be combined with the equations from Eqs. 
([5]) - ([7]) to obtain the full system of equations that governs 
the phase dynamics of the diamond- shaped bi-SQUID 
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Fig. 4. Calculated average voltage response of a single diamond shaped DC 
SQUID, performed with i b = 2.0, L\ = 0.54, L 2a = L 2b = L 5a = L 5b = 
0.24, L 3a = L 3b = L 6a = L 6b = 0.3 and i c3 = 0.5. 




Fig. 5. Diamond- shaped dual bi-SQUID: (a) cell schematic; (b) micropho- 
tograph of the fabricated cell within a 2D array. I c = 0.25 mA, = 2.4 
fl , V c = I c R s h = 600 fiY. All junctions are critically shunted, j3 c = 1. 



Ln{<Pi ~ 03) + ^(206 + 05 - 04) 



Li 



^26 ib 



+ ^2 - <Pi + 27rx e ai + L12 (i C 3 sin <p 3 - sin (^1) 
+ (sin (f 4 - sin y? 5 - 2i c6 sin <p 6 ) 



Li 



(20 6 + 5 - 4 ) - L 12 ((f2 + 03 ) 



2a fc& 



+ y? 2 - y>i + 27rx e ai + L12 (i C 3 sin <p 3 + sin <p 2 ) 
+ ^ (sin <p 4 - sin <p 5 - 2z c6 sin <p 6 ) 
^303 = ^2 - </?i - ^3 - £3^3 sin (p 3 

^15(04 - 06) + ^(203 + 02 - 0l) = (^Lhb + ^6 

+ y?5 — y?4 + 27nr e a 2 + L15 (i C 6 sin y? 6 - sin <£> 4 ) 



+ ^ (sin (px - sin <p 2 - 2z c3 sin <p 3 ) 



Li 



(20 3 + 02 - 0l) - ^15(05 + 06) 



^56 + -y J ^ 



•^5 - ^4 + 27rx e a 2 + L15 (i c6 sin <p 6 + sin <p 5 ) 



+ y (sin - sin <p 2 - 2z c3 sin (p 3 ) 
L 6 6 = y> 5 - (p 4 - cp 6 - L 6 i c3 sin (p 3 , 



(8) 



where L 15 = L\ + ^5a + £56, ^6 = ^6a + ^66 and we use 
the approximate assumption that a 2 = L\ + L 5a + L55. 

The calculated average voltage response for the diamond 
shape is shown in Fig. [4] This average voltage response is 
more linear than the single SQUID but less so than the single 
bi-SQUID, shown in Fig. [2] The non-linearities are most likely 
due to the inclusion of additional inductances in the diamond 
structure of bi-SQUIDs. The analysis is now extended to the 
full 2D diamond array. The phase equations that model the 
2D diamond arrays are derived in a similar way as the single 
diamond structure of bi-SQUIDs. These modeling equations 
total 66. Both the system of equations and the derivations are 
not expressed explicitly for brevity. 

IV. Implementation of 2D Diamond bi-SQUID 
Arrays 

For experimental evaluation of our 2D SQIF arrays of 
diamond- shaped dual bi-SQUID cells, various size arrays were 



laid out. Fig. [5] shows a schematic and a layout of the diamond- 
shaped dual bi-SQUID cell designed for a standard HYPRES 
Nb fabrication process |43]-|45). R s h is the shunting resis- 
tance and V c is the critical voltage across each of the junctions. 
Its electrical circuit was extracted and re- simulated to account 
for the actual layout parameters. The layout was done using 
all four available Nb layers: a ground plane layer, two layers 
for junctions and inductors, and a top layer to implement a 
flux bias line overlaying bi-SQUID cells (two light gray strips 
retracing the diamond cell contour). The ground plane was 
used only under Nb layers forming bi-SQUID inductors and 
junctions in order to maintain their low specific inductance. 
However, ground plane was partially removed from under 
inductor L\ to increase its value if desired. The ground plane 
was also removed from the central area of the bi-SQUID loops 
to allow an external magnetic field to thread through the cell 
(dark areas in Fig. [5jb)). The total area of the dual-bi- SQUID 
cell is ~ 1623 /inr . While the total area of the cell is kept 
the same, the values of inductances are varied by their width 
and ground plane opening size under L\. We designed arrays 
with normal Gaussian distribution of cell inductances with a 
~ 30% and ~ 70% for comparison. 

Fig. [6] shows a schematic and the layout implementation 
for a 2D array using diamond dual bi-SQUID cells. It is 
evident that a diamond shape of our dual bi-SQUID cells leads 
to their natural 2D arrangement into a diamond checkered 
pattern by connecting dual bi-SQUID cells by their corners. 
The space between cells forms a similar array comprising 
dual bi-SQUIDs with somewhat different ratio of inductors 
forming bi-SQUID loops. As a result, each bi-SQUID shares 
its junctions (Ji and J 2) and inductances with neighboring 
bi-SQUIDs. Each dual bi-SQUID has a contact to three 
neighboring cells at each of four corners (8 in total). This 2D 
array design avoids the use of long parasitic wires, as every 
component of the array is an essential element of a bi-SQUID 
indicating to an efficient use of the available area. 

The array dc bias is fed uniformly from the top bi-SQUID 
row of the array. Similarly, the array is grounded to the 
bottom bi-SQUID row. The inductively coupled flux bias line 
is overlayed on the top of the array forming loops for each 
column. The direction of the dc flux bias control current is 
shown by small arrows in Fig. |6jb). This line can be used for 
RF signal input for testing in the low noise amplifier regime, 
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while the output signal is measured at the dc bias current 
terminal via a bias tee. 




Fig. 6. 2D diamond- shaped bi- SQUID array: 
microphotograph of the fabricated 2D array. 



(a) 2D array schematic; (b) 



Fig. [7] shows microphotographs of some examples of the 
fabricated 2D bi- SQUID SQIF arrays with different size and 
configuration. We chose to resistively shunt all junctions in 
the bi-SQUIDs in the arrays shown. The objective of these 
designs is to investigate the array dimension dependences. 
Only 12 contact pads are used. These are concentrated on 
one side of the chip in preparation for antenna tests in which 
RF signal will be irradiated to the chip. To investigate the 
array characteristics, several test arrays with different sizes 
and configurations were laid out using 5x5 mm 2 chips for 
the fabrication with a 4.5 kA/cm 2 Josephson junction critical 
current density. 

V. Experimental Evaluation as a Low Noise 
Amplifier 

Experimental evaluation was first performed in liquid he- 
lium using the test setup described in detail in [46]. Each 
chip was tested using HYPRES standard cryoprobes with 40 
coaxial lines. These cryoprobes can make contact only to the 
required smaller number (12) of pads (Fig. [7} without shorting 
the array structure by remaining 38 unused cryoprobe contacts. 

A. Experimental Investigation of Noise Properties 

Fig. [8ja) shows the measured flux-to- voltage characteristics 
of the 15 x 80 bi- SQUID array with a 70% inductance spread 
shown in Fig. [7ja). As one can see, it has a well-defined 
zero-field anti-peak. Fig. [8jb) shows the corresponding flux 
noise measurement at the mid-point of the positive slope of 
the anti-peak. The flux noise spectral density is reaching to ~ 
2 x 10 _6< l ) o/V / Hz, which is the expected value for this array. 
Fig. [9] shows the corresponding energy sensitivity and noise 
temperature calculated following (47), (48). 

The noise energy per unit bandwidth via flux noise in a 
SQUID is 

«*(/) 





(a) (b) (c) 

Fig. 7. Microphotographs of 5 x 5 mm 2 chips of different 2D diamond bi- 
SQUID SQIF arrays with a a ~ 70% spread: (a) a 1200 bi-SQUID (15 x 80) 
array; (b) two serially connected 2D arrays (2 x 43 x 85) arrays with 7310 
bi-SQUIDs; (c) a single 7820 bi-SQUID (92 x 85) 2D array. The sizes of 
these arrays are calculated by assuming each diamond has two bi-SQUIDs. 
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Diamond BiSQUID Array (Chip#12-M352-1) 




(b) 



10 J 10' 
Frequency (Hz) 



Fig. 8. Measured characteristics of a 15 x 80 cell dual bi-SQUID SQIF array 
of Fig. [7J a) with a ~ 70% of inductance spread: (a) flux/voltage characteristic 
with 2 mV/div 0.5 mA/div, max voltage w 18 mV, AV/AI (flux bias) ^170 
V/A. (b) measured flux noise spectral density. The spikes are attributed to a 
noisy frequency generator. 



where / is frequency. The inductance, L, of bi-SQUID is 
calculated from the measured separately AI C modulations of 
the IV curve defined as L = . The noise temperature can 
be expressed as 

*fe(f) 



(10) 



e(f) 



2L 



(9) 



where ks is Boltzmann's constant. The quantum limit of 
noise temperature is defined as Tql = \hf /ks |49|. The 
upper frequency (100 kHz) for our sensitivity measurements 
was determined by the available test equipment. However, the 
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Diamond Bi SQUID Array (Chip#12-M352-1) 




(a) 



10 10 10' 10 J 

Diamond BiSQUID Array (Chip#12-M352-1) 



(b) 




10" 10 

Frequency (Hz) 



Fig. 9. Noise characteristics of the 1200 (15 x 80) bi-SQUID SQIF array 
shown in Fig. [7|a): (a) energy sensitivity; (b) noise temperature compared 
against quantum limit. 



noise temperature has a linear frequency dependence for given 
device, e. g., for a SQUID RF amplifier, T/v scales as the 
ratio ujq/V(<&) (50). So, one can extrapolate the noise data 
obtained at lower frequencies (100 kHz) to higher frequencies 
and estimate the noise properties up to, e.g., ~ 2 GHz, which 
is still much lower than Josephson frequency at the DC bias 
point. 

B. Experimental Results Compared Against Simulations 

For comparison to results of simulation described above, the 
arrays chips were fabricated with a 30% and 70% Gaussian 



spread in inductances. Figs. [TO] and [TT] show that increasing the 
Gaussian spread to 70% makes the average voltage response 
go from a wave packet-like response to a single anti-peak. The 
single anti-peak is the form of the average voltage response 
in which we are interested. The generation of the single anti- 
peak feature is significant because the anti-peak was obtained 
with identical size bi-SQUID loops and with varying the 
inductances. The fabricated arrays that produced the current 



flux/voltage characteristic in Figs. 10 and 11 were 15 x 40 
and 15 x 80 arrays, respectively. While the simulations were 
performed with arrays of size 3 x 80. To simulate larger arrays 
the code will need to be parallelized. 

A comparison of the experimental and simulated current 
flux/voltage characteristic was performed on an array of size 




(a) 

Fig. 10. Comparison of the (a) measured (5 mV/div, 10 mA/div) for a 
15 x 40 dual bi-SQUID array and (b) simulated flux (control current)/voltage 
characteristic for a ~ 30%. The simulations were performed using a 3 x 80 
array with = 2.0 and i C 3 = 1.0. 




Fig. 11. Comparison of the (a) measured (2 mV/div, 1 mA/div) for a 
15 x 80 dual bi-SQUID array and (b) simulated flux (control current)/voltage 
characteristic for a ~ 70%. The simulations were performed using a 3 x 80 
array with = 2.0 and i C 3 = 1.0. 




I $ (mA) 

Fig. 12. Comparison of the measured and simulated control current 
flux/voltage characteristic for a 15 x 20 diamond array with the inductances 
spread a ~ 30%. The measurements of the fabricated array were performed 
with Ibias = 2.4 and 2.8. The simulations were performed with = 1.5 
and 1.75 and i C 3 = 0.5. 



15 x 20, see Fig. 12 While there are many similarities, the 
difference in the width of the peaks is most likely due to the 
array having unshunted third Josephson junctions while the 
system of equations that model the system are derived using 
all resistively shunted junctions. 

C. Experimental LNA response 

Measurements of the input power and output power for sin- 
gle tone frequencies were performed 5.3 MHz, 208 MHz, and 
413 MHz at the CERF Center of Excellence on a 7 x 80 array 
with a 30% Gaussian spread in inductances and unshunted 
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TABLE I 

SQIF LNA Output to Single RF CW Tones 



Input Power 5.3 Mhz (dBm) 


Normalized Power Output (dB) 


-100 


_ 


-80 


25 


-60 


42.5 


Input Power 208 Mhz (dBm) 


Normalized Power Output (dB) 


-100 





-80 


14 


-60 


34 


Input Power 413 Mhz (dBm) 


Normalized Power Output (dB) 


-100 




-80 


15 


-60 


34.5 



Table of power input and output values to a single tone at frequencies 5.3 
MHz, 208 MHz, and 413 MHz. 
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Fig. 13. Visualization of the points of device bias and flux bias for 
measurements. 



third junctions, see Table |T| The measurements were performed 
with the output power normalized to the lowest observable 
power, which was about —126.4 dBm, corresponding to an 
input power of —100 dBm at 208 MHz, current bias of 1.3 
mA, and flux bias from a control line current of 0.54 mA, 



see Fig. 13 More measurements will need to be performed to 
determine how the device is performing as an LNA, however 
the data will be useful when characterizing devices in the near 
future (for RF studies). 

D. Antenna Sensitivity Analysis 

We estimated possible antenna sensitivities for a SQIF- 
based antenna by assuming that the area for the diamond 
shaped (double) bi-SQUID is 1.62 x 10 -9 m 2 . Using that as 
the effective area, multiplying by 15 x 80/2 to account for the 
number of individual bi-SQUIDs in the array, and taking the 
flux noise from Fig. ^b) as 2 x 10 -6 &o/VRz, gives a field 
sensitivity of ~ 4.25 fT/v^Hz. 

Assuming a scaling as y/N, the field noise for a 1000 x 2000 
array would be 0.104 fT/y/Hz at 100 kHz. Approximating the 
physical dimension of the diamond as ~ 71 fim x 71 /im 
means that the diamond occupies an area of ~ 5 x 10 -9 m 2 



(consistent with the effective area of the diamond shaped bi- 
SQUID). Thus a 1000 x 2000 array would occupy an area of 
~ 50 cm 2 , since each diamond contains two bi-SQUIDs. This 
corresponds to a square with 7.1 cm on a side. 

VI. Conclusion 

We have introduced a new design of 2D SQIF arrays using 
tightly-coupled diamond- shaped dual bi-SQUID cells. This 
design allowed us to address the challenge of preserving 
and enhancing the linearity of individual bi-SQUIDs which 
tends to be degraded when bi-SQUID cells are connected in 
parallel. Another advantage of this approach is that the cells 
can be densely packed onto a chip, resulting in many more 
bi-SQUIDs in the same area when compared with strings of 
bi-SQUIDs coupled in series. Our SQIF array design involves 
bi-SQUIDs with identical areas and varying inductances. We 
were able to show that a anti-peak feature can be generated 
with these features. This approach may be applicable for 
arraying conventional SQUIDs as well. 

A system of equations was derived to model the dynamics 
of the 2D diamond dual bi-SQUID array. The model was very 
complex and to model large arrays the Matlab code will need 
to be parallelized. We found a good match between simulations 
of small arrays and experimental results. The array design 
has undergone extensive computer simulations, analysis and 
experimental measurements, in which we explored the voltage 
response of the 2D arrays. 

Experimental measurements were performed on the noise 
properties of the array design. The results matched the ex- 
pected values. More extensive experimentation with larger 
size arrays of different configurations are required in order to 
optimize the array design to meet the requirements of wide 
band low noise amplifiers and electrically small antennas. 
Ongoing work includes Transverse ElectroMagnetic (TEM) 
cell measurements to experimentally evaluate the array design 
as an antenna. 
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